% sobel算子
img=imread("img/liftingbody.png");
subplot(3,3,2),imshow(img),xlabel("原始图像");

MaskPrewittV=fspecial("sobel");
MaskPrewittH=MaskPrewittV';

KB1=filter2(MaskPrewittH,img);
subplot(3,3,4),imshow(uint8(KB1)),xlabel("水平模板滤波图像");
subplot(3,3,5),imshow(uint8(double(img)+KB1)),xlabel("水平模板滤波加法叠加图像");
subplot(3,3,6),imshow(uint8(double(img)-KB1)),xlabel("水平模板滤波减法叠加图像");

KB2=filter2(MaskPrewittV,img);
subplot(3,3,7),imshow(uint8(KB2)),xlabel("垂直模板滤波图像");
subplot(3,3,8),imshow(uint8(double(img)+KB2)),xlabel("垂直模板滤波加法叠加图像");
subplot(3,3,9),imshow(uint8(double(img)-KB2)),xlabel("垂直模板滤波减法叠加图像");